Antiferromagnetic to valence-bond-soild transitions in two-dimensional SU(N) 
Heisenberg models with multi-spin interactions 
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We study two-dimensional Heisenberg antiferromagnets with additional multi-spin interactions 
which can drive the system into a valence-bond solid state. For standard SU(2) spins, we consider 
both four- and six-spin interactions. We find continuous quantum phase transitions with the same 
critical exponents. Extending the symmetry to SU(iV), we also find continuous transitions for iV = 3 
and 4. In addition, we also study quantitatively the cross-over of the order-parameter symmetry 
from Z4 deep inside the valence-bond-solid phase to U(l) as the phase transition is approached. 

PACS numbers: 75.10.Jm, 75.10.Nr, 75.40.Mg, 75.40.Cx 



Two-dimensional quantum spin system with non- 
magnetic ground states have been at the forefront of 
condensed matter physics for more than two decades 
[2 0, H J3| • Frustrated system have been investigated in- 
tensly 5| , but large-scale unbiased computational studies 
of their ground states are not possible, due to the "sign 
problems" hampering quantum Monte Carlo (QMC) 
methods [fjj. It was recently realized that one promi- 
nent class of non-magnetic states — valence-bond solids 
(VBSs) — can be accessed also without frustration, by 
adding certain multi-spin interactions to the standard 
S = 1/2 Heisenberg antiferromagnet 0- These mod- 
els enable detailed QMC studies of the antiferromagnetic 
(AF) to VBS quantum phase transition. It has been 
argued that this transition is associated with spinon de- 
confinement (hence the term deconfined quantum criti- 
cality) and should, due to subtle quantum interference ef- 
fects, be continuous 0]. This scenario violates the "Lan- 
dau rule" , according to which a direct transition between 
states breaking unrelated symmetries should be generi- 
cally first-order. 




FIG. 1: (Color online) Interactions involving p singlet projec- 
tors (illustrated by ovals enclosing two sites) on the square 
lattice. The two-spin (p = 1) interaction J is the Heisenberg 
exchange. Higher-order Q p terms with p = 2 and 3 are con- 
sidered here. All translations and 90° rotations of the site 
groupings shown here are included in the hamiltonian. 



is defined by the hamiltonian 

Hi = J Si ■ Sj = — J Cij 

(y) <*J> 



L 2 J 



(1) 



where (ij) denotes nearest neighbors on a periodic square 
lattice with L 2 sites and 



1 



Si • Sj 



(2) 



The theory of deconfined quantum criticality has gen- 
erated a great deal of interest^ as well as controversy 
0, S, 0, EE EH, 0, 0, El- Numerical studies of 



is the two-spin singlet projector. In the "J-Q" model 
introduced in Q, the following term is added to Hi; 



a Heisenberg hamiltonian with 4-spin interactions are 
generally in good agreement with the theory, showing 
a continuous transition with dynamic exponent z = 1, 
large spin correlation exponent rj s , and an emergent U(l) 
symmetry @, @, • Arguments for a first-order transition 
have also been put forward ll|, [L4| , based on numerical 
studies of lattice versions of the field-theory proposed [3( 
to capture the AF-VBS transition. Other, similar studies 
reach different conclusions, however (l3| . Further studies 
are thus called for. 

In this Letter, we advance computational studies of the 
AF-VBS transition in two different ways. First, we con- 
sider the S = 1/2 Heisenberg model including 4-spin and 
6-spin interactions. The unperturbed Heisenberg model 



H% — —Q2 ^ C'kiCij. 

(ijki) 



(3) 



The spin pairs ij and kl are located on adjacent corners of 
a 4-site plaquette, as illustrated in Fig.[TJ We denote the 
strength of the 4-spin term Q2, with the subscript indi- 
cating two singlet projectors, and also consider a similar 
term with three stacked singlet projectors, 



H3 — —Q3 ^ CmnCklCij , 
(ijklmn) 



(4) 



as also illustrated in Fig. [TJ Using an improved version 
[IE of a ground-state QMC method operating in the 
valence-bond basis 17[, we have studied the J-Q2 and 
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FIG. 2: (Color online) Finite-size scaling of the squared AF 
and VBS order parameters of the J-Q2 and J-Q3 models. 

J-Q3 models on lattices with L up to 64. We find crit- 
ical AF-VBS points with the same set of exponents for 
both models, providing additional evidence of a universal 
deconfined critical point in this class of systems. 

In a second development, we have studied SU(N) sym- 
metric versions of the J-Q2 model, in the representation 
of the spin operators previously used in mean-field 0] and 
QMC calculations Ij3 of the SU(N) Heisenberg model. 
We find continuous AF-VBS transitions also for N = 3 
and 4 (whereas for N > 4 the system is VBS ordered for 
all Q 2 > [HEll). 

An open problem in previous studies of the J-Q2 model 
was that the order parameter distribution inside the VBS 
phase did not show the expected 4-fold symmetry. In- 
stead, the distribution was always U(l) symmetric 0, 
An emergent U(l) symmetry close to criticality is in- 
deed predicted by the field theory Q as a consequence 
of a dangerously irrelevant operator, but deep inside the 
VBS phase the order parameter should exhibit Z4 sym- 
metry (which has been observed in other quantum models 
0,0). With the J-Q 3 model and the N > 2 versions 
of the J-Q2 model, we can now reach sufficiently deep 
inside the VBS phase to observe the expected U(l)-Z4 
cross-over. We present quantitative finite-size scaling re- 
sults for the exponent governing the cross-over. 

For all the models, we compute the square of the stag- 
gered magnetization, M 2 = (M • M) , where 

1 



M = ;p£(-ir 



,r,y 



V &x,y (5) 

is the operator of the AF (spin) order parameter. We 



model, symmetry 


r) s 


Vd 


V 


0,4 


J-Q2, 
J-Q3, 


SU{2) 
SU{2) 


0.35(2) 
0.33(2) 


0.20(2) 
0.20(2) 


0.67(1) 
0.69(2) 


1.20(5) 


J-Qa, 


517(3) 


0.38(3) 


0.42(3) 


0.65(3) 


1.6(2) 


J-Qa, 


517(4) 


0.42(5) 


0.64(5) 


0.70(2) 


1.5(2) 



TABLE I: Critical exponent for all the models studied. The 
cross-over exponent 04 cannot be determined for the SU(2) 
J-Q2 model, because no cross-over is observed for L < 64. 

define the columnar VBS order parameter in terms of 
neartest-neighbor (dimer) correlators 



(6) 



and D y defined analogously. We compute the square 
D 2 = (D 2 + Dy) and also study the probability distribu- 
tion P{D X , D y ), with D x and D y evaluated in the config- 
urations generated in the QMC sampling (as in 0]). To 
extract the critical points and exponents, we use stan- 
dard finite-size scaling forms for the order parameters, 



M 2 = L 1 +"'F,([g- 9c ]L 1 /") i 
D 2 = L L+Vd Fd([q — q^L 1 ^), 



(7) 

(8) 



where ry s and r\d are the exponents governing the spin 
and dimer correlation functions, respectively, at critical- 
ity (the anomalous dimensions) and 1 + r\ s _d = 2/3 s ,d/ 1 '. 
Here we assume a dynamic exponent z = 1, in accord 
with previous studies of the J-Q2 model @, S] , and use a 
single correlation length exponent v, as in the theory 0] • 

We first present results for the SU(2) models. Defining 
coupling ratios q = Q V /{J + Q P ), we find critical points 
q c = 0.961(1) for p = 2 and q c = 0.600(5) for p = 3. The 
former agrees with previous estimates @, S, 0] • Standard 
data collapse plots according to ([7|) and (|HJ) are shown 
in Fig. O The critical exponents are listed on the first 
two lines of Table HI Here it is very significant that all 
the exponents are the same for the two models. This 
supports the notion of a universal deconfined quantum 
critical point. Note that the order parameters decay as 
L 1+,,s ci at the common critical point q = q c . At a first- 
order transition, the order parameters should instead be 
size-independent at q c , due to phase coexistence. 

Comparing with previous results for the J-Q2 model, 
the results for smaller systems in 0] were consistent with 
Vs = Vd (with a value between those found here), but 
the present results for larger systems clearly show that 
the spin and dimer exponents are different. The theory 
does not make any specific predictions for a relationship 
between r) s and rjd, and they can be expected to be dif- 
ferent. The exponents rj s and v are in good agreement 
with values obtained using finite-temperature scaling in 
[8( (where ^ was not determined). 

Next, we discuss the J-Q2 model generalized to SU(JV) 
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FIG. 3: (Color online) Scaling of the spin and dimer order 
parameters of the SU(3) and SU(4) J-Q2 models. 

spins. Considering first the Heisenberg model, the hamil- 
tonian can be written as 



H 



SU(7V) 



N 



0a 
3 



(ij) 



2JL 2 
N 2 



(9) 



where S" 13 is the generator of the SU(iV) algebra, with 
a, (3 — 1,2, - ■ ■ ,N the different "colors", and is the 
generalization of Eq. ([2|) to S\J(N). As in we focus 
on the simplest case, where the spins on sublattice A are 
expressed in the fundamental representation (i.e., with 
a single-box Young tableau). Spins on sublattice B are 
SU(iV) conjugates (dual representation) of those on A (a 
Young tableau with one column and N — 1 rows). The 
states in this representation can be written in terms of 
permutations P of the boxes, with 



—EC- 1 ) 



'\P(2)P(3)---P(N)) j , 



(10) 

with a = 1, 2, • • • , N and P(l) = a. An SU(A) singlet of 
spins i and j on different sublattices is given by 



I singlet) i 



1 



N 

Q — 1 



(11) 



QMC algorithms using these SU(A) spins in the valence- 
bond basis are simple generalizations of the SU(2) case 
0, E3, EH- Instead of spins | and j. for SU(2), there 
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are N colors, and, thus, N states of the space-time loops 



FIG. 4: (Color online) Dimer order distribution P(D X , D y ) 
for L = 32 systems. The left panels are for the J-Q3 model 
at q = 0.635 (a) and q = 0.85 (b), and the right panels are 
for the SU(3) J-Q 2 model at q = 0.45 (c) and q = 0.65 (d). 



in the loop-algorithm [16|. The off-diagonal matrix ele- 
ments of the singlet projection operators are 1/N instead 
of 1/2, and the overlap of two valence-bond states is gen- 
eralized to N n °~ L / 2 , where n is the number of loops 
in the transposition graph. The 4- and 6-spin terms ([3]) 
and ([4]) are written explicitly using products of singlet 
projectors and have obvious generalizations to SU(iV). 

Our results for the SU(3) and SU(4) versions of the J- 
Q2 model are consistent with continuous AF-VBS crit- 
ical points, with no signs of first-order behavior. The 
critical points are q c = 0.335(2) and q c — 0.082(2) for 
N = 3 and 4, respectively. Scaling plots giving the crit- 
ical exponents are shown in Fig. [3] and numerical values 
are listed in Tablefl] It can be noted that, as a function of 
N, v does not change appreciably, rj s increase slowly, and 
rjd increases significantly. An increasing spin exponent is 
consistent with rj s = 1 in the N — * 00 theory As 
already noted, there are no predictions for rjd- We could, 
in principle, consider still higher A, but with J > the 
system is always in the VBS state for N — 5 and higher 
[lil [l^ | . A transition could presumably be reached for 
J < 0, but this causes QMC sign problems. 

The dimer order distribution P{D x ,D y ) can be used 
to investigate the VBS order parameter symmetry 
As shown in Fig. [H for large q the robust VBSs in the 
SU(2) J-Q 3 model and the SU(3) and SU(4) versions of 
the J-Q2 model result in histograms with clearly visible 
columnar Z4 features (i.e., peaks on the D x and D y axis, 
as opposed to 45° degree rotated histograms expected for 
a plaquette state). However, in the SU(2) J-Q2 model the 
histograms are ring-shaped for all system sizes currently 
accessible, even in the extreme case of q = 1 (J = 0). 
In all cases, we see U(l) symmetric histograms as the 
critical point is approached, in agreement with one of 
the salient features of deconfined quantum criticality 0] . 
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FIG. 5: (Color online) Finite-size scaling of the square of the 
anisotropic order parameter D4 in the J-Q3 model (upper 
panel) and SU(3) J-Q2 model (lower panel). 

Defining an order parameter sensitive to the symmetry, 

D\ = J dD x dD y P(D x ,D y )(D 2 x +D 2 y )cos(4e) 



dr I d8P(r,6)r A cos(4(9), 




(12) 



where 9 is the an gle corresponding to a point (D x ,D y ), 
we proceed as in [22j (which deals with a classical sys- 
tem with a dangerously irrelevant perturbation) to ex- 
tract the exponent governing the length-scale A of the 
Z4-U(l) cross-over (and the spinon confinement). Z4 fea- 
tures should appear for L > A, which is predicated Q 
to scale as A ~ £ a4 where £ is the correlation length and 
a 4 > 1. We analyze D4 assuming the scaling form [2^ ] ; 



Dl = L 1 - 



F A (qL l/aiV ). 



(13) 



This form describes the cross-over, as shown in Fig. [5] in 
two cases. The values of 04 are listed in Table HI The 
large error bars reflect slow evolution of the VBS angle 
in the QMC simulations. It is nevertheless clear that 
04 > 1 (and increasing with TV), reflecting emergent U(l) 
symmetry due to a dangerously irrelevant perturbation. 

The results presented here support deconfined quan- 
tum criticality. Although one can still, in principle, not 
completely rule out very weakly first-order transitions 
based on these calculations, the universal behavior for 
the two SU(2) models makes this less likely. The com- 
mon exponents for the J-Q2 and J-Q3 models at the very 
least suggest close proximity to a universal critical point. 



The detailed information now available from QMC sim- 
ulations should be useful to further advance the theory. 

In a very interesting experimental development, Itou 
et al. have recently measured the spin-lattice relaxation 
rate 1/Ti in a la yer ed organic compound which seems to 
be near-critical [23j. It has been argued that, in spite 
of the triangular lattice, the AF-VBS transition in this 
kind of system should be in the same class of deconfined 
quantum critical points discussed here 24j. The expo- 
nent rj s governs the temperature scaling of 1/Ti, and the 
value rj s « 0.35 is in excellent agreement with the exper- 
iment over a wide range of temperatures. Further exper- 
iments should elucidate the nature of the ground state 
and whether it is indeed close to a deconfined quantum- 
critical point. 
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